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Abstract. A short introduction to the theory of matrix quasiclassical Green's func- 
tions is given and possible applications of this theory to transport properties of meso- 
scopic superconducting- normal metal (S/N) structures are considered. We discuss a 
simplified version of these equations in the diffusive regime and in the case of a weak 
proximity effect. These equations are used for the calculation of the conductance of 
different S/N structures and for analysis of kinetic phenomena in these structures. We 
discuss the subgap conductance measured in SIN tunnel junctions and the mechanism 
of a nonmonotonic dependence of the conductance of a N wire on temperature T and 
voltage V, observed in an S/N structure. 

Long-range, phase- coherent effects are studied in a 4-terminal S/N/S structure un- 
der conditions when the Josephson critical current is negligible (the distance between 
superconductors is much larger then the coherence length in the normal wire). It is 
shown that the Josephson effects may be observed in this system if a current /, in 
addition to a current I\ in the S/N/S circuit, flows through the N electrode. 



1 . Introduction 

The progress of nanotechnology over the last few years has made it possible to 
produce conducting nanostructures in which new physical phenomena have been 
observed. Specifically, hybrid structures consisting of superconductors (S) and nor- 
mal conductors (N) have been created using metal films [2-5] or semiconductor 
layers [1,6-8] as the normal conductors. The transport properties of these S/N 
structures have turned out to be quite unusual. First, the subgap conductance 
(zero-bias anomaly) has been observed in SIN tunnel junctions at low tempera- 
tures (T <100 mK) [l](see also [7,8]). Second, conductivity oscillations have been 
observed in these mesoscopic structures in a magnetic field H (i.e., in structures 
with dimensions less then the phase-breaking length L^). Oscillations of the con- 



ductivity of the N channel appeared if the structure contained superconducting or 
normal loops [2-6]. Moreover, for an N channel in contact with superconductors a 
nonmonotonic dependence of the conductivity on the temperature T and voltage V 
has been observed at T^CT C [5]. The main experimental facts have been explained 
in recent theoretical works (see review articles [23,24]). It was established that the 
proximity effect plays the main role in the transport properties. For example, the 
conductivity of an N channel in the structure shown in Fig.l changes as a result 
of the contribution of the condensate induced by the proximity effect. Since the 
condensate is induced by both superconductors in a nonlocal manner, interference 
appears and a term — 5R cos ip, which depends on the phase difference f between 
the superconductors, arises in the resistance of the N channel [9-11]. The phase 
difference increases with the magnetic field H, and this results in oscillations of the 
conductivity of the N channel in a magnetic field. The nonmonotonic dependence 
of the resistance R of an N channel on T and V has also been explained [12-15] 
(see also the theoretical works in the Conference Proceedings in Ref. [14]). The 
nonmonotonic dependence of the resistance R(T, V) of a point contact ScN (c is a 
constriction) was first obtained theoretically in Ref. [16]. 

New effects have also been predicted in theoretical work devoted to S/N struc- 
tures. For example, in Refs. [14,17] it was shown that the critical Josephson current 
I c in a structure of the type displayed in Fig.l depends on the voltage Vg between 
the S and N conductors, changing sign ( n - contact) if Vs exceeds a certain value. 
In addition, it has been shown that the Josephson effect also arises in the case when 
current flows only through one S/N boundary. Several different configurations of 
S/N structures were studied in Ref. [18], determining that under certain conditions 
the current- voltage characteristics of the S/N structures have descending segments 
(dl/dV < 0). 

An important circumstance was noted in Ref. [19] (see also the works in Ref. [14]). 
It was shown that the local conductivity of an N channel changes over distances 
from the S/N boundary which can be much greater than the coherence length 

£at = \J D/2itT (D is the diffusion coefficient). Important consequences follow from 
this fact. For example, phase coherence effects in the conductivity of an N channel 
remain even if the distance 2L\ between the superconductors is much greater than 
£at. This means that the conductivity oscillations in the structure shown in Fig. lb 
will also be observed in the case of a negligibly low critical current I c . The oscillation 
conservation effect is due to fact that as T increases, I c decreases exponentially 
(I c ~ exp(-2Li/^ iV CO), and 5R decreases slowly (5R ~ T~ l ) [20,35]. 

In these lectures we discuss briefly the method of quasiclassical Green's functions 
and apply this method to the study of transport phenomena in mesoscopic S/N 
structures. We restrict ourselves to the dirty limit where the mean free path / is 
essentially less than geometric dimensions of the system and the coherence length, 
but exceeds significantly the Fermi wave length hp (quasiclassic approximation). 
We will consider mostly a weak proximity effect, when the amplitude of the conden- 
sate induced in the normal metal is small compared to the condensate amplitude 



in the superconductors S. This case occurs if the S/N interface resistance is larger 
than the resistance of the normal conductor N. Results obtained for this case remain 
qualitatively valid in case when these resistances are comparable. 

In the next Section we present the main necessary equations for the Green's func- 
tions and a general expression for the current in the N channel in which a condensate 
is present due to proximity effect. In Section 3 we will give formulas describing the 
subgap conductance of tunnel S/I/N junctions and discuss a possible physical in- 
terpretation of this conductance. In Section 4 we will consider the conductance of 
a N channel attached to two superconductors and obtain a formula describing, in 
particular, the oscillatory behaviour of the conductance in an applied magnetic field 
H . Also a nonmonotonic dependence of the conductance on temperature T and on 
bias voltage V will be analysed. The possibility of observing Josephson-like effects 
in a S/N/S mesoscopic structure (see Fig. lb) will be considered in Section 5. We 
will show that zero voltage between superconductors may exist in some interval of 
the current through the S/N interfaces and Shapiro-steps may be observed even in 
absence of the real Josephson coupling between superconductors when the distance 



a) 



-V 




V 


N 




1 

S. 




N 


N I 




s + 
1 





~L 



w w b) 







1 s 


Il II 


1 s 


















s. 




s + 






N 












N 


I 


N I+l! 




I 




A 





Li ~t 



FIGURE 1. Schematic diagram of the system considered 



separating superconductors exceeds essentially the coherence length 



2L 1 >f JV (T) = ^D/2ttT. (1) 

It is important that these effects arise only in the case when a current / flows along 
the N channel and the dissipation takes place [21,35]. 

2. Basic equations for quasiclassical Green's functions 
The Green's function technique is a powerful tool for the theoretical study of differ- 
ent phenomena in superconductors and superconducting structures. In the case of 
superconducting systems, we need to indroduce condensate Green's functions of the 
type < (1) ipi (2) >, therefore all the Green's functions have a matrix form. For 
example, the retarded (advanced) Green's functions are defined as follows [22,23] 

G R W = ±9(t m - t 2{1) ) (G> (1,2)- G< (1, 2)) (2) 
Here G K and CP* are 

G> fi = -i< &,(1WJ(2) > (-1) Q+1 , 

G< /3 = ,<^ + (2)^(l)>(-ir +1 (3) 

We introduced here spin indices in the Nambu-space: ipi (1) = ^ (1), tp 2 (1) = 
(1). As is well known, the functions G R ^ describe the excitation spectrum of 
the system. In order to describe nonequilibrium processes, one needs to know the 
distribution functions which are related to the Green's function G indtroduced by 
Keldysh. It is convenient to define a supermatrix G (1,2) elements of which are the 
matrices G R ^ = (6(1,2)) and G = (6) . The element (g) is the zero 

V V ' ' ) 11,(22) V /12 V /21 

matrix. 

In the quasiclassical approximation all components of the Green's functions 
G ex (1, 2) are integrated over the variable £ p = {p — Pf)vf and in the corresponding 
equation for G ex (1, 2) an expansion is carried out in the parameters (pird) _1 ,(pirZ) _1 
or (pf^n)" 1 , where d is the thickness of the S or N films, £ N is the coherence length 
in the N conductor and / is the mean free path. The quasiclassical Green's functions 
are defined by the relation [22,23] 

G (p/p, f; t u t 2 ) = (i/n) J d£ p G ex (p, f; t ± , t 2 ) (4) 

The subscript "ex" means exact (nonquasiclassical) Green's functions. Therefore 
the quasiclassical Green's functions G (1, 2) depends on the angle of momentum on 
the Fermi surface, on the coordinate r, and on two times. 

In what follows we need an equation for the supermatrix (5(1,2) only in the N 
conductor, having the form 

DV (GVG) + ie [a z , 6\ = 0. (5) 



where a z is a diagonal supermatrix elements of which are the Pauli matrix a z . 
Eq. (5) may be averaged over the thickness of the N film d. Performing the averag- 
ing, we obtain 



Dd x (Gd x G) + le [a z , 6\ = e b 9{x s ) [G s , G . 



(6) 



where the coefficient e& is a characteristic energy which is proportional to the trans- 
mission of the S/N boundary: e& = pD/2Rbnd, Rm is the resistance of a unit area 
of the S/N boundary; p and d are the resistivity and thickness of the N film. 
When deriving Eq. (6), the boundary condition 



was used. Here the z-axis is normal to the plane of the S/N interface. The bound- 
ary conditions for the quasiclassical Green's functions G have been derived in the 
general case by Zaitsev [25] and have been reduced to the simple form (7) by 
Kupriyanov and Lukichev [26] in the dirty case. In the case of a good S/N con- 
tact, condition (7) is reduced to the continuity of the Green's functions at the S/N 
interface: G = Gs- In the case of a poor contact (e& — > 0), condition (7) gives the 
same result for the current through the S/N interface as obtained with the aid of 
the tunneling Hamiltonian method. However, for a S/N contact with an arbitrary 
barrier transparency condition (7) is not applicable. The point is that when de- 
riving Eq. (7) Kupriyanov and Lukichev [26] restricted themselves to the Legendre 
polynomials of the zeroth and first orders in the expansion of the angle-dependent 
Green's function G. Meanwhile, one can easily show that all the Legendre har- 
monics are excited near the S/N (or N/N') interface. They decay to zero (except 
the Legandre polynomials of the zeroth and first order) over the mean free path 
away from the interface. In order to obtain a correction of the next order in e& to 
condition (7), one has to solve an integral equation (see Ref. [27]). In the case of 
the S/N interface with an arbitrary barrier transparency, the problem of boundary 
conditions for the quasiclassical Green's functions becomes complicate (boundary 
conditions and their applicability are discussed in detail in the Raimondi's Lecture 
Notes). 

Eq.(6) must be solved in the normal conductor for a particular geometry (see, 
for example, Fig.l) with boundary conditions at x — ±L. In the case of normal 
reservoirs the condensate functions F R ^ (±L) are equal to zero and G R<yA) (±L) = 
±a z . In the case of superconducting reservoirs the boundary conditions for the 
retarded (advanced) Green's functions are 



where G Ri - A) = e/^ A \ = A/^ R{A) (±ia x sin<f> + ia y cos<j>) ,^ {A) = 



(e±ir) 2 — A 2 , T is a damping in the excitation spectrum in superconductor, 



20 is the phase difference between superconductors. Eq.(8) is valid if the voltage 




(7) 



8 




between superconductors 2V is much less than A/e. We assume that there is no 
barrier between the N conductor and reservoirs. 

The Keldysh function G describes the kinetic properties of the system and is 
related to distribution functions 

G = G R f - fG A (9) 

where / = f Q l + fa z is a matrix distribution function. The function f a enters 
an expression for the supercurrent (in a superconductor it determines the energy 
gap), and the function / determines the quasiparticle current (in a superconductor 
it describes the charge- imbalance and the electric field; see, for example [28]). In 
reservoirs the functions /„ and / are supposed to have equilibrium form 

fo(±L) = [tanh((e + eV)(3) + tanh((e - eV)(3)]/2, (10) 
f(±L) = ±F N (e) = ±[tanh((e + eV)(3) - tanh((e - eV)(3)]/2, (11) 

where f3 = {2T)-\ 

If the functions G R ^ and G are known, one can easily find a relation between 
the applied voltage 2V and the current / in the N conductor. The expression for 
the current is [22,23] 

/ = {ad/8)Tra z J de(G R d x G + Gd x G A ) (12) 

Eqs.(6) and (12) can be simplified significantly in the case of a weak proximity 
effect when the amplitudes of the condensate functions in the N conductor F R ^ 
are small. Then the retarded (advanced) Green's functions in the N conductor have 
the form 

QR(A) =G R(A) dz+ p R{ A) (13) 

where G R{A) w ±[1 + ( y F R( - A ^) 2 /2] and | F R W |< 1. When obtaining the relation 
between G R ^ and F R ^ A \ we employed the normalization condition [22,23] 

(G R W) 2 1 + (F R ^) 2 = 1, (14) 

where 1 is a unit matrix. 

The equation for the condensate functions F R ^ follows from Eq.(6) and the 
expression (14) 

d xx F R{A) - (k R ^) 2 F R ^ = -k 2 b w {P*l A) 5 (x - L ± ) + F^ A) 5 (x + L,)] (15) 

where k\ = 2e b /D, (k Ri - A ^ 2 = (±2ie + i)/D, w is the width of the S/N interface in 
the x-direction and is supposed to be much less than L and £jv(T). We introduce 
here the depairing rate 7 in the N conductor which determines the phase-breaking 



length L§ = y.D/7. Once the functions F R ^ are known from a solution of the 
linear equation (15), we can find the conductance of the N film. 

Let us consider for example the system shown in Fig.l. Writing down Eq.(6) for 
the matrix element (12), we arrive at the equation outside the S/N interface 

Dd x [& R d x G + Gd x G A ] + ie [a z , 6] = 0. (16) 

Multiplying Eq.(16) by a z and calculating the trace, we obtain after the first inte- 
gration 

(0 x /)[l-m_] = J(e). (17) 
Here J(e) is an ^-independent constant and 

m- = (\/$)Tr(F R -F A f (18) 

is a function which describes the condensate contribution to the N film conductance. 
The left side of Eq.(17) stems from the first term in the square brackets in Eq.(16) 
provided that the condensate functions are small. Therefore, according to Eq.(12), 
the current I is an integral from the "partial current" J 



{ad/2) J deJ(e) (19) 



Solving Eq.(17) with boundary conditions (11), we can find a relationship between 
the current and voltage I(V). In the next Sections we will analyse the conductance 
of S/N mesoscopic systems. 

3. Subgap conductance in SIN junctions 
In this Section the subgap conductance in superconductor/insulator/normal metal 
(S/I/N) tunnel junctions will be discussed. As is well known from conventional 
theory for S/I/N junctions, the subgap conductance should exponentially decrease 
with decreasing temperature T (see, for example, Ref. [29]). However, experiments 
on Nb/n+InGaAs contacts have established that a peak in conductance appears at 
zero-bias if the temperature becomes low enough (T < A), and the magnitude of 
this peak at low temperatures (T ps 50 mK) is comparable with the conductance 
in the normal state [1]. This contact can be considered as a tunnel S/I/N junction. 
A Schottky barrier at the interface plays the role of the insulating layer I. An 
explanation for anomalous transparency of the SIN junction at low voltages and 
temperatures (T,eV <C A) was suggested in Refs. [9,30-34]. According to the 
interpretation proposed in Ref. [32], the subgap conductance is due to a component 
of the current which, in the case of a SIS Josephson junction, gives the so-called 
interference current. This component can be presented in the form 

I mt = - (8R,)- 1 j de-F N (e, V) (f r + F A ) (f r + F A ) , (20) 



where F N (e, V) is defined in Eq. (11); F* (A) = A/ ((e ± iY) 2 - A 2 ) 1/2 are the 



condensate, retarded (advanced) Green's functions in the superconductor. This 
formula can be obtained from the general expression for the current (12) and from 
Eq.(6) in which we have to put 6{xs) = 1- If the energy e is small, F R = F A m 
—i. In the case of a S/I/N junction, F R W 

are the condensate functions in the N 
electrode. To zero order in the barrier transmittance (i.e., in -R^ 1 ), they are equal 
to zero. If the proximity effect is taken into account, they differ from zero and in 
the case of a planar S/I/N junction they have the form (see, for example, Ref. [32]) 



F R(A) 



±e b / (e ± ry) , 7 > e b 

1/2 21 



e b / (£±17) -e 



This formula can be obtained from Eq.(6) in the case of weak and strong proximity 
effect. It is seen from Eq.(21) that F R ^ are small if e b <C £,7, where e is deter- 
mined either by temperature T or voltage V . In the opposite limit when e and 7 
are small compared with e b , F R ^ are not small, and the differential conductance 
normalized by R^ 1 , S = R b dI/dV, calculated from Eq.(20) for T=0 and V=0 is 
not small either. The integrand in Eq.(20) is not zero if |e| < A because F R = F A 
at |e| < A and F R = —F A at \e\ > A. This means that the current given by Eq.(l) 
is caused by the charge-transfer mechanism of the same type as Andreev reflection 
processes. The second important circumstance leading to the subgap conductance 
is related to an anomalous proximity effect when the amplitude of the condensate 
functions F R ^ at small energies e is not small. 

The density-of-states (DOS) in the N electrode is changed drastically in the 
case of the strong proximity effect: the DOS is zero at | e |< e b and the DOS = 

e/^/e 2 — e b in the interval A e |> e b . In a one-dimensional S/I/N junction the 
DOS has a quasigap at | e |< e b . In both cases of a planar or one-dimensional S/I/N 
junctions the zero-bias, zero-temperature conductance coincides with it's value in 
the normal state [32]. 

4- Conductance of the Andreev interferometer 
Consider the system shown in Fig.l (the Andreev interferometer). In order to 
calculate the normalized differential conductance of the N channel S = R^dl/vV in 
the presence of a phase difference between superconductors, we must solve Eq.(17) 
taking into account the boundary condition (11). The function m_ is small by 
assumption. Therefore the expression for J may be presented in the form 

J(e)=F N (e,V)[l-(m-)]/L (22) 

where (m_) = (1/L) J L dx ■ m 2 _ . Substituting (22) into Eq.(19), we can obtain 
an expression for S. Here we present the formula for a deviation of the nor- 
malized differential conductance from it's value in the normal state R^: 5S = 
(2L/ad) dl/dV - 1 . We obtain 



5S = - (1/2) Jde-(3-F' N (e, V) (m_) , 



(23) 



where F N (e, V) = [cosh" 2 (e + eV) (3 + cosh' 2 (e - eV) (3 
inition of (mJ) (see Eqs.(18)), we can write (m_) in the 



/2. By virtue of the def- 



orm 



(m_) = Tr((P R ) 2 + (P A ) 2 - 2F R F A j /8 (24) 

The first two terms in (24) determine a change in the DOS of the N channel due 
to the condensate (this term reduces the conductance), and the last, so-called 
anomalous, term leads to an increase of the conductance. 

As it is seen from Eq.(23), in order to find the conductance, we need to solve 
Eq.(15) for the condensate functions F R ^ (x). In this Section we present here the 
solution for the geometry shown in Fig. la. 

F R W{x) ^ia y F^ A) r[e R ^coshe R ^Y\mh[k R ^(L-\x\)]-cos V (25) 

Here r = k^Lw is the ratio of the N channel resistance to the S/N resistance, 
9 R(A) = k R(A) L _ Calculating (m_) we find 

(m_) = (V 2 /8) [Re [sinh (26) /26 - 1] / (0cosh0) 2 - [sinh (29 Y ) /26 1 

- sin (20 2 ) /20 2 ] / \0 coth 6\ 2 } (1 + cos 2ip) (26) 

where 6 = Q\ +i9 2 - The first term in (26) determines a contribution to the conduc- 
tance due to a change in the DOS, and the second term is related to the anomalous 
term (F A F R ). In Fig.2 we show the dependence SSjjos (first term contribution) 
and the 5S an (V) dependence (anomalous term contribution) at T = 0. It is seen 
that 5Sdos (V) is negative and SS an (V) is positive. The total change in the con- 
ductance 5S (V) = SSdos iV) + 8S an (V) is shown by the solid line. This quantity 
increases with increasing V from zero, reaches a maximum at V m ~ e^/e , and 
decays to zero with further increase of V (here ei = D/L 2 is the so-called Thouless 
energy) . 

As follows from (23) and (26), the conductance SS oscillates with increasing the 
phase difference. In Fig.2 we also plot the temperature dependence of the zero-bias 
conductance. Both curves, 5S (V) at T = and SS (T) at V = are similar. 
Note that nonmonotonous temperature dependence of the conductance was ob- 
tained earlier in Ref. [16] where a short ScN contact was analyzed (here c means a 
constriction). 

5.Dissipative Josephson-like effects in S/N/S structures 
In this Section we discuss a possibility to observe Josephson-like effects in meso- 
scopic S/N/S structures (see Fig. lb) with negligible Josephson coupling between 
superconductors, i.e., when the inequality (1) is fulfilled [21,35]. Following the same 
steps as in Section 2, we obtain instead of Eq.(17) 



(1 -m_ (x))d x f 



J + Ji - J 5 , < x < Li 
J, Li < x < L 



(27) 



In what follows the function m_ plays the most important role. 

The current on the segment (0, L±) is determined by Eq.(19) if J is replaced 
by J + J\. The quantity J s , the superconducting "current", is constant over the 
segment (L±, L) and (0, L\) and is equal to 



J s = (1/4) Tra z (P R d x F R - F A d x F J 



(2f 



The integral of Js (28) over the energy is exponentially small if the condition (1) 
is satisfied. As follows from Eq. (6), the constant J\, is related with the Green's 
function and distribution function in the superconductor. It can be written in the 
form [11] 



Here 



Ji 

R bD /w 



J q + Js, J q = [p/d N U b ) [F s (e) - / {L x )\ 



(29) 

is the resistance of 



~v N v s + (1/8) Tr [P R + F A ) (Pf + F t 
the S/N boundary per unit length in the y direction and un, vs - are the density 
of states in the N and S conductors. It can be shown that for Vn,s which are small 
compared with T/e, the "supercurrent" Jg, flowing throw the S/N boundary equals 
Js- The distribution function Fs is the equilibrium function, i.e., it is identical to 
the function in Eq. (11), if Vn is replaced by Vs (we measure voltages from the 
point 0, where the voltage is equal to zero). Using the fact that m_ is small, we 
can integrate Eq. (27) and find the relation of J and J q with F^ and Fs (see the 
boundary condition (11)). We obtain the normal currents 



(d N /p) J = 
(d N /p) Ji w J q (d N /p) 



M b F N + (F N - F s ) 
dt 2 Fs + 3?i (F s — F N ) 



(30) 




FIGURE 2. Normalized conductance SS vs. normalized voltage eV/e^ at T = and vs. nor- 
malized temperature T/cl at V = for the structure shown in Fig. la. 



Here % is determined in Eq. (29); the quantity 3? = 3?i+9? 2 , ^1,2 = Ri,2 (1 + ( m -)) 
can be termed the partial resistance. The spatial average (m_)i j2 on the segments 
(0, Li) and (L 1 , L) gives a decrease in the resistances on account of proximity 
effect ((rri-) is negative). All resistances in Eq. (30) depend on the difference of the 
phases ip and on the energy; they can be represented in the form 9fy, = Rb — 5$t,b cos ip 
and 3?i )2 = -Ri,2 — $fti,2 cosy?. The corrections to the resistances <53f? 6 and <55ft 1)2 are 
small in the case of a weak proximity effect. The quantities Rb and Ri^, depend, 
generally speaking, on the energy e (for example, v$ depends on e). We assume, 
for simplicity, that these quantities do not depend on the energy. This is valid if it 
is assumed that the superconductors are gapless (the results remain qualitatively 
the same in the case of superconductors with a gap). Then, integrating Eq. (30) 
over energies, we obtain on the left-hand side the currents I and l\ (see Eq. (19)). 
Eliminating Vn from the two equations obtained, we find for Vs 

V s = hd t ip/2e = 

h [R b + Ri- {5R b + SRx) cos 2ip] + I (R 1 - SRi cos 2ip) (31) 

Here we employed the Josephson relation; Rb is the resistance of the S/N boundary, 
which in the case of zero-gap superconductors is approximately equal to its value in 
the normal state. The resistance R\ is also approximately equal to pLi/d,N (the f 
-independent correction arising from (m_) is small and unimportant). Integrating 
Eq. (31), we obtain a relation between the average voltage Vs and the fixed currents 
I and 1 1. 

V s = y/[(I + h) Ri + hRbf - [(/ + h) SRi + hSR b } 2 (32) 

The function Vs (h) is displayed in Fig. 3 for different currents /. One can see that 
for I 7^ this dependence is identical to the current-voltage characteristic of a 
standard Josephson contact. In this case the critical current is 

T _ T 5R\Rb — 5RbR\ , . 

Therefore I c increases in proportion to the current I. We shall show below that 
the correction 5R\ decreases slowly with increasing temperature {5R\ ~ T" 1 ), and 
the correction 5R b is small if the condition (1) is satisfied. Therefore, for R b ^>R\, 
we obtain I c ~ ISRi/Rb- The maximum current / is limited by the condition that 
Joule heating be small and by the condition eV/v — elR^T . In the opposite case 
5R\ decreases as Vjy increases. If the condition (1) is not satisfied and a finite 
Josephson coupling exists between the superconductors, then it is easy to show 



that the critical current of the structure equals I* = + where I c j is the 
critical Josephson current. An expression for I c j can be easily obtained with the 
aid of Eq. (28). This expression is presented in Ref. [20]. The equilibrium phase 
difference ipo for Ii + lR\j [Rb + R\) = equals 2ip = — arcsin(/ c //*). 



To determine 5Ri and SRb it is necessary to find the condensate functions F R ( A \ 
For \x\ < L\ the solution of Eq.(15) has the form 



F R(A) (x) = F R(A) [ia y cos (<p) P y cosh (kx) +ia x sin (<p) P x sinh(fca;)] 



R{A) 



(34) 



R(A) 



■ s is the amplitude of the condensate functions in the superconductors. 
In the zero-gap case Fg^ = ±A/ (e ± ijs), where 75 is the frequency of spin-flip 
collisions with impurities. The functions P XyV equal: 
P x = b sinh 62 1 (sinh 6 + b sinh 9\ sinh 9 2 ), 
P y = b sinh 62/ (cosh 9 + b cosh 9\ sinh 9 2 ), 

b = pw/ (R bn d N ) k, k R ^ = y/T2ie/D, 9 = 9 1 + 9 2 , 61,2 = 0[ >2 + %9'[ 2 = kL 1>2 Once 

the functions F R ^ are known, the interference correction 5R\ to the resistance 
can be calculated: 

6R1 = -RiJ de(3 ■ cosh _2 (e/3)(m_ (x,(p) - m_ [x,n/2)) 1 (35) 


With the aid of the expressions for (m._)i and for F R ( A ' (see Eq. (34)), we find 

00 

5Ri/R x = J de(3- cosh- 2 (e/5)M (e) , (36) 



where M (e) = (1/8) [\F S 

- \P X \ 2 [sinh [29[) /29[ - sin (261) 



\P y \ 2 [sinh (20[) /29[ + sin (29'l) /20?] 
+ ReF§ [P 2 (sinh (200 /2^ + 1) 




FIGURE 3. Vs versus the current I\ for the following values of the current: 1 - 0, 2 - 250/iA, 3 
- 500/uA, 4 - 750^A, 5 - 1mA. Here 5R X = O.LRi, i? 6 = 5i?i, i?i = IO. 



— P% (sinh (26 1) /29i — 1)] j. The temperature dependence of 5R\ is displayed in 

Fig. 4. One can see that for T > = D/ (2Li) 2 the quantity 5R\ decreases as 
T _1 with increasing temperature. As noted in Refs. [15,18], the slow decrease of 
SRi (T) is due to the so-called anomalous term F r Fa in (m_)x. The special role of 
this term, which is nonanalytic both in the upper and lower planes of e, was noted 
in Ref. [36]. 

The Josephson current Is is determined by the integral of J,s (28), over all 
energies, i.e., the integral of products of either advanced or retarded Green's 
functions. It can be calculated by closing the integration contour in the upper 
(lower) half plane of e and switching to summation over the Matsubara frequencies 
uj n = 7rT (2n + 1). For such energies the functions F R ^ decay exponentially over 
distances k~ x (u n ) < £ n (T) away from the S/N boundary. Therefore the current 
Is will be exponentially small (Is ~ exp (— 2Li/^ (T)). The function Is{T) for 
the structure shown in Fig. lb is presented in Ref. [20]. Similar arguments are also 
applicable to the calculation of SRf,, since for T < 75 the functions Fg and Fg 
can be assumed to be equal and independent of the energy. At the same time, 
the function F R F A , appearing in the expression for 8R1, decreases over a small 
(compared with T) energy = Dj (2L X ) 2 and makes a nonzero contribution. For 
such energies the characteristic decay length of F R ( A >(x) is of the order L 1 , i.e., of 
the order of the distance between the superconductors. 

In order to observe long-range Josephson effects, the critical current I c must 
exceed the fluctuation current Te/h: I c ~^>Te/h. On the other hand, the ordinary 
Josephson effect is negligible if the condition e^^cT is fulfilled. Combining these 
inequalities, we obtain the condition 

TR b R 1 /{5R b R Q )<.e I ^<r (37) 




FIGURE 4. Interference correction 8R\ to the resistance as a function of temperature in the 
case L x = Q.5L, R/R b = 0.4, 7 /e L = 100, A/e L = 30. 



which should be satisfied for observation of the effects under consideration. Here 
Rq = h/e 2 « 3/cfi, and we took into account that a maximal value of / is deter- 
mined by the relation eIR < e^. Otherwise SRi decreases with increasing /. The 
first inequality of (37) means that the zeroth Shapiro step on the I\ (V5) curve is 
absent at / = 0. If the second inequality of (37) is not fulfilled, then the critical 
current is not zero at / = (ordinary Josephson effect). In this case the effective 
critical current I* should first increase with increasing / and then decrease when I 
exceeds e^/eR. 

6. Conclusion 

In conclusion we note that, as one can see from Fig. 4, the correction 8R\ to the 
resistance of the normal channel caused by the proximity effect depends on the 
temperature T in a nonmonotonic way: it is equal to zero at T = (the bias 
voltage is zero as well), reaches a maximum at T and decays to zero at higher 

T. Such behavior of 8R1 (T) is related, as noted in [15], to different dependencies 
of two contributions to 5Ri on the energy e. One contribution which increases 
the N channel resistance is connected with a decrease of the density-of-states in the 
normal channel, which is described by the last term in M (e) (see Eq. (36)). Another 
contribution (anomalous) which diminishes the resistance of the normal channel is 
described by the first two terms in M (e). This contribution exactly compensates a 
contribution due to a change in the density-of-states of the normal channel at e = 
and dominates at e 7^ 0. At T > T c it leads to the Maki-Thompson contribution 
to the paraconductivity. Mathematically, compensation of the two contributions 
at e = arises because at e = F R = F A and m_ in Eq. (35) tends to zero. The 
nonmonotonic behavior of 5R has been observed in an experiment [5]. It would be 
interesting to observe the long-range Josephson effect experimentally. 
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